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Abstract 

Bohmian mechanics, a hydrodynamic formulation of quantum mechanics, relies on the concept of 
trajectory, which evolves in time in compliance with dynamical information conveyed by the wave 
function. Here this appealing idea is considered to analyze both classical and quantum interference, 
thus providing an alternative and more intuitive framework to understand the time-evolution of waves, 
either in terms of the flow of energy (for mechanical waves, sound waves, electromagnetic waves, for 
instance) or, analogously, the flow of probability (quantum waves), respectively. Furthermore, this 
procedure also supplies a more robust explanation of interference phenomena, which currently is only 
based on the superposition principle. That is, while this principle only describes how different waves 
combine and what effects these combinations may lead to, flow lines provide a more precise explanation 
on how the energy or probability propagate in space before, during and after the combination of such 
waves, without dealing with them separately (i.e., the combination or superposition is taken as a 
whole). In this sense, concepts such as constructive and destructive interference, typically associated 
with the superposition principle, physically correspond to more or less dense swarms of (energy or 
probability) flow lines, respectively. A direct consequence of this description is that, when considering 
the distribution of electromagnetic energy flow lines behind two slits, each one covered by a differently 
oriented polarizer, it is naturally found that external observers’ information on the slit crossed by single 
photons (understood as energy parcels) is totally irrelevant for the existence of interference fringes, in 
striking contrast with what is commonly stated and taught. 

Keywords: Bohmian mechanics, flow line, interference, sound wave, electromagnetic wave, Umov vector, 
Poynting vector, probability current 

1 Introduction 

Superposition and interference are two intertwined capital concepts in both classical and quantum 
wave mechanics. As is well known, within the framework of a general wave theory, superposition simply 
consists in adding two or more waves. In elementary courses on wave mechanics [lj|^, it is commonly 
taught that when two traveling waves propagate through a certain medium, the net displacement of the 
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medium at a particular point and a given time is just the direct sum of the individual displacements 
associated with each contributing wave. 

Such a traditional conception of interference, where only the net effect of the waves is considered 
and the associated flow of energy is often neglected (or just pushed into the background), strongly relies 
on the mathematical grounds of wave mechanics, namely the superposition principle. According to this 
principle, given a set of field-like variables {i9i}, solutions to a linear flow equation 

= 0 , ( 1 ) 

where denotes the (space and time) flow operator 

^ ^2 2 ^2 

=-c — 

dx^ ’ 

and c is a constant typically identified with the propagation (or diffusion) 
solutions can be readily built up from linear combinations of the 

N 

Q = '^Ci'di. ( 3 ) 

i=l 

This intuitive principle has an important practical advantage: it allows us to find and devise smart 
methods to solve linear flow equations, like Eq. Q, by just decomposing the total solution into a set of 
partial contributions. Besides, in our opinion it also brings in an important inconvenience: our conception 
of physical superposition is directly grounded on its mathematical basis, which becomes critical when 
Eq. 0 represents a wave equation and the field variables refer to waves propagating from different 
sources. As a consequence, for example, following the above mathematical roots, it is commonly thought 
that when two waves superimpose, they pass each other just as if each wave was unaware of the other; only 
the net effect matters [^. This is just a physical model to understand and explain the wave phenomena 
observed in Nature. However, it also strongly determines our perception of such phenomena, what has 
very important consequences particularly in the case of quantum systems. 

To some extent, such a picture of interference is rather poor. It refers to the superposition or combi¬ 
nation of partial waves, although in real life these waves never appear separately. In order to get some 
more light on interference, let us get back to the usual notion of wave from classical wave theory, where 
it typically represents a perturbation that propagates throughout a certain medium. Erom a physical 
viewpoint, the propagation of a perturbation in a material medium is associated with a flow of energy 
from one place to another within such a medium. A better understanding of interference should there¬ 
fore include the tools to determine how this flow takes place, and not only rely on how different wave 
components (which we never observe in real life) combine. The same idea, of course, can be naturally 
extended to quantum systems. Even if the nature of quantum waves is different from that of classical 
waves, the key element, namely the transport or flow of a certain quantity (probability in the case of the 
former and energy for the latter) is the same. 

In this work we tackle the issue of interference within a general theoretical framework based on the 
concept of flow line, applicable to both classical and quantum systems. Within this scenario, while the 
superposition principle only tells us how different waves combine, the trajectories or flow lines provide 
us with information on how the net effect propagates before, during and after their combination, putting 
the emphasis on the more natural conception of considering the wave as a whole, and not on splitting it 
up in its different components, as it is typically done. 


( 2 ) 


speed of the field variable, new 
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It is worth noting that, although the idea of transversal flow is mentioned in different sources (for 
example, see energy flow for non-dispersive waves in [^, or for the probability flow in quantum processes 
in 0[|)- the idea of monitoring it during the full propagation of the wave (i.e., time by time) is not that 
general at all. As far as we know, the use of flow lines as a working tool to visualize the propagation of 
energy in sound waves can be traced back to the 1980s (6 - 10 . The concept of flow line is quite general and 


can be found in many different physical contexts as well [11| , including classical electromagnetic problems 
(see Ref. 12 and references therein). Specifically, in most of these works this concept is introduced after 


assuming an analogy between the corresponding problem and a hydrodynamic one. 

Together with the fact that flow lines constitute an interesting working tool to analyze the evolution 
(flow) of energy or probabilities, we would also like to stress that depending on the problem considered, 
these elements present the nice feature that they can be somehow inferred from experimental data, as 
shown recently in the case of Young’s two-slit experiment 
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This representation is in compliance 
with the recent method devised to determine the photon wave function [14] , based on the so-called 
weak-measurement technique [^, which constitutes a remarkable alternative to the more traditional 
construction of phase-space tomograms [16] . 


This paper is organized as follows. In Sec. we discuss sound waves as an example of classical waves, 
introducing the associated equation for the energy flow lines. As an illustration, we present a simulation 
of the energy flow lines behind a wall with two openings. In this example, the nodes in the eventual 
fringe pattern represent regions of silence (or, to be more precise, low sound) in between regions where 
the sound is reinforced. This experience could be equally reproduced with two synchronized loudspeakers 
emitting the same sound pitch. The homologous equation for electromagnetic energy (EME) flow lines 
is given in Sec. Moreover, in analogy with the previous case, in this section we also show the EME 
flow lines behind two slits, although each slit is followed by a linear polarizer with its polarization axis 
oriented in a different direction. In Sec. |^we discuss how an even more general picture of interference 
can be achieved for quantum nonzero-mass particles after considering flow lines of quantum-mechanical 
probability density. Einally, the main conclusions arising from this work are summarized in Sec. 


2 Sound Energy Flow Lines: Interference behind Two Openings 


In order to derive the equation of motion accounting for the transport of energy of a sound wave propa¬ 
gating through a non-viscous fluid, we start from the continuity and momentum equations [3,17,18 , 


dp 


+ V (pv) = 0, 


5v 

p ——h Vp + pv • Vv = 0. 
ot 


(4) 

(5) 


respectively, where p(r, t) is the fluid density, v(r, t) is the velocity of a fluid element, and p(r, t) is the 
pressure. Assuming the velocity of the fluid and the fluctuations of the density and pressure are small 
(p = Po + P^ p' ^ PO) P = Po +P^ p' ^ PO) V = v'), after linearization the above equations become 


|-+pVv' = 0. 


( 6 ) 

(7) 
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Taking into account the relationship between pressure and density, 


dp 

dp 


= ^-. = c\ 


(8) 


P=PO 


where c is the speed of sound in the corresponding medium, from Eqs. and 0 we readily obtain the 
wave equation for the acoustic pressure. 


The energy transport in the sound field is described by the equation 


m 


T»'’'' + 2^)+^k'v')=o, 


(9) 


( 10 ) 


where 


S = p'v' (11) 

is the so-called Umov energy flow vector [3,19 , the mechanical equivalent of the Poynting or Poynting- 
Heaviside vector, commonly used in electromagnetism [20] . 

Let us now consider without loss of generality and for simplicity a sound wave of constant angular 
frequency uj (the analysis could be carried out as well with non-monochromatic waves). If the pressure and 
velocity are taken as complex quantities (a typical working technique when dealing with wave equations, 
although only the real part of these quantities is physically meaningful). 


p' = 

after substitution into Eqs. 0 and Q, we have 


v = Ve 


—iujt 


v^p+ 0) P = 0, 


iwV = 


VP 
Po ' 


( 12 ) 

(13) 

(14) 


In the case of high-frequency waves, to avoid fast oscillations, we can consider time-averaging over a 
cycle, which allows us to introduce the definition of the time-averaged energy flow vector. 


{S) = ^Re{PV*} 


The modulus of this vector determines the sound intensity |17| 

i = m\, 

related to the energy density w as 

c c 


(15) 


(16) 


(17) 


Consider now the propagation of the above monochromatic plane sound wave incident onto an obstacle 
with two slits, located on the XZ plane. Assuming the amplitude of the incident complex pressure does 
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Figure 1: Energy flow lines associated with sound waves emerging from an obstacle with two slits. The 
width of the openings is 5 = 0.5 m, their mutual distance d = 1.5 m, the speed of sound c = 340 m/s, 
and its frequency / = 1 kHz. 


not depend on the z coordinate and has the functional form the solution of the Helmholtz 


equation (13) behind the obstacle can be expressed in terms of the well-known Fresnel-Klrchhoff integral 


from the theory of propagation of optical waves 20 - 22 


P(x,y) = r 

V /_oo 


(18) 


Here, P{x\ 0) is the initial pressure, which in our case is given by the pressure just behind the obstacle. 


determined by the corresponding boundary conditions. From Eqs. (14) and (15), we thus obtain 


1 f dP* dP 

(S) = - 7 ^- Im <! P ^ e,, + P ^ e 


2poUJ 


dx 


dy 


y ( > 


(19) 


where ex and ey are unit vectors along the x and y directions. The vector field (S) is tangent to the 
sound energy flow lines at each point determined by the equation 


dr I 

Behind a one-dimensional grating, this equation can be recast as 

dx {Sx)' 


( 20 ) 


( 21 ) 


Sound energy flow lines calculated from Eqs. (18)-(21) for a totally absorptive obstacle with two 


identical openings are shown in Fig. We find that the energy flow lines associated with different slits 
do not cross the system symmetry line. This means that there is a kind of fictitious barrier between the 
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two sound sources (slits). The flux-line non-crossing property for classical/sound waves is analogous to 
the Bohmian non-crossing property pS-25 . This correspondence can be set on the basis of the analogy 


between Eq. ( |20[ ) and the equation of trajectories determined by the quantum probability current density 
(see Sec. [^for a more detailed discussion). 

It is worth noticing that sound waves are described by field variables, i.e., continuous functions of 
space and time, although the gases, liquids or solids through which these waves propagate consists of 
individual, discrete objects (e.g., atoms, ions, molecules, etc.). However, in the case of fluids, for instance, 
instead of appealing to these objects to analyze the corresponding hydrodynamics, it is common to make 
use of tracer particles, which help to elucidate the flow dynamics and hence the dynamical properties of 
the fluid 


26,27 


The propagation of sound waves (energy) is of practical interest, particularly in room acoustics [^. In 
this area, sound energy flow lines have been evaluated in the past using different approximate methods. 
Recently, the diffraction of these waves at half planes and slits has been studied p9 30 by means of 


methods of scalar optics in the far field. The method described above also makes use of results from 
optics (in general, electromagnetism), although is more general, since it encompasses both near and far 
field. 


3 Electromagnetic Energy Flow Lines: Interference behind Two Slits 
Covered by Polarizers 

Can the electromagnetic field associated with light be thought as a kind of “wave function for the photon” ? 
Although this is a longstanding challenging question, it is still of much interest nowadays. Scully and 
Zubairy summarized 31 the objections posed by Kramers |32], Power 33 and Bohm 34 against the 


concept of a photon, agreeing with Kramers and Bohm that the concept of a photon wave function should 
be carefully used, since it can be very misleading. Nonetheless, they concluded that each objection can 
be overcome, suggesting a path to circumvent them, and providing a way to construct the photon wave 
function through the radiation-field second quantization formulation. 

Raymer and Smith 


35 , Bialinicki-Birula 36-^ and Holland [39| went a step further. Holland 
analyzed similarity of Maxwells’ equation with quantum wave equations (Schrodinger’s and Dirac’s) 
Raymer and Smith 


35 and Bialinicki-Birula 136 -38 


consider that the usual electromagnetic Maxwell 
field is the quantum wave function of the single photon. The fact that it transforms like a three- 
dimensional vector arises from the spin-one nature of the photon. The interpretation of the Maxwell 
field is, therefore, akin to the Schrdinger wave function, which describes the evolution of probability 
amplitudes associated with various possible quantum events in which a particle’s position is found to be 
within a certain volume. In this sense, Maxwell’s equations rule the evolution of probability amplitudes 
for various possible quantum events in which the photon’s energy is found within a certain volume. 


Assuming this interpretation for the Maxwell field, and based on a previous work by Prosser 140,41 


a method to determine the electromagnetic energy flow lines behind various gratings illuminated by a 


monochromaticbeam of light was developed in 12,42 . In this case, the solutions of Maxwell’s equations 


behind the gratings were determined using the solution of the Helmholtz equation. (It is well known that 
the space-dependent part of each component of the electric E(r) and magnetic H(r) free fields satisfy 
Helmholtz’s equation). 

As before, again here we are going to consider the time-averaged energy flux vector to determine the 
electromagnetic energy flow lines (for a non-averaged, time-dependent application, see the work by Chou 
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and Wyatt 1^), given by the real part of the complex Poynting vector 44 


S(r) = Re 


-E(r)xH*(r) 


( 22 ) 


Since the energy flow is along the direction of the Poynting vector, the EME flow lines will be determined 
from the parametric differential equation 

dr S(r) 

ds cC/(r)’ 

where s is a certain arc-length and U (r) is the time-averaged electromagnetic energy density 


(23) 


U{r) = - [eoE(r) • E*(r) + /ioH(r) • H*(r)] . 


(24) 


As mentioned in the previous section, the Poynting vector in electromagnetism plays the same role as 
the Umov vector (19) in the theory of mechanical waves. Since in the last instance both vectors are 


just transport vectors, some authors use the combined denomination Poynting-Umov vector (or also 
Heaviside-Poynting-Umov vector, to include Heaviside’s contribution). 

Holland argued that, for reasons of compatibility with quantum mechanics, it does not seem 
to be reasonable assuming that the tracks r = r(t) deduced from the equation essentially equivalent to 
Eq. (23) can be identified with the orbits of’’photons”. However, Ghose et al. showed |45| that Bohmian 


trajectories for relativistic bosons, and so for photons, could be defined indeed. They derived an equation 
of motion for massless bosons, which is equivalent to Eq. (23). From such an equation they obtained 


photon trajectories behind the two-slit grating, assuming that the wave function at the point {y,x), at 
the sufficient distance D ^ d'^/X to the right of the plane of the slits, is a superposition of two spherical 
waves. 

EME flow lines evaluated from Eq. (23) in |12[[^ show that the energy redistribution behind a multiple 
slit grating corresponds with a Talbot pattern in the near field, and with a Fraunhofer interference pattern 
in the far field. The Fresnel-Arago laws as well as the Poisson-Arago spot were also interpreted using 
EME flow lines 46,47 . The main conclusion from the analysis of these phenomena is that the motion 


of an eventual photon wave packet essentially represents the flow of electromagnetic energy along an 
ensemble of flow lines. 

Actually, EME flow lines obtained |47[|48| from a numerical simulation of Young’s two-slit experiment 
with parameters taken from the experiment performed by Kocsis et al. [I^ showed a good agreement 
with the averaged photon trajectories inferred from the experimental data. It is worth stressing that the 
photon trajectories reconstructed from the experiment were in compliance with the Bohmian approach, 
thus confirming that trajectories coming from different slits do not cross. This means that, at the level of 
the average electromagnetic field (or the wave function, in the case of material particles, in general), full 
which-way information can still be inferred without destroying the interference pattern. That is, rather 
than complementarity, the experiment seems to suggest that a photon wave function has a tangible 
(measurable) physical reality [^, in agreement with a recent theorem on the realistic nature of the wave 
function [49| . 

Let us now evaluate the EME flow lines behind two slits in a more general case, when each slit is 
covered by a polarizer with its polarization axis oriented at a different angle. In our model, we consider a 
monochromatic electromagnetic wave in vacuum, traveling along the y direction, incident onto a two-slit 
grating located on the XZ plane, at y = 0, with infinitely long slits parallel to the 2 ; axis, so that both 
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the electric and magnetic fields do not depend on the ^ coordinate. The electric and magnetic fields of 
the incident wave are given by the expressions 


E, 


Hi, 




^J^o 




(25) 

(26) 


Since the slits are followed by polarizers, with their polarization axis at angles 9i and 62 with respect to 
the z axis, the ^ components of the fields behind the grating read as 


E2(r) = ^cos^ 0iV^i(r) + Tcos^ 02V’2(r) 

—Be^y^ sin 01 cos0iV’i(r) — sin 02 cos02V’2(r) 


H2(r) = —Acos 0 i sin 0 i'i/’i(r) — ^cos02 sin02V’2(r) 


eo 






+ \l — Be^y’ sin^ 0iV'i(r) + \ l — Be'’'^ sin^ 02V'2(r). 


Mo 


Mo 


(27) 


(28) 


Where V’i(i') s-iid V' 2 (r) are scalar functions that satisfy the Helmholtz equation and the boundary con¬ 
ditions at the slits. The functions V’i('') V’ 2 (r) can be recast in the form of a Fresnel-Kirchhoff 

integral. 






(29) 


with i = 1,2, and where the '^/’^(^c^0^) denote their value just behind the first and the second slits, 
respectively, and in front of the polarizers. 

Since the electric and magnetic fields do not depend on the z coordinate, from Maxwell’s equations 
one obtains two independent sets of equations j^: one involving the Hx and Hy components of the 
magnetic held and the component of the electric held (commonly referred to as E-polarization), 


dE, 

dy 


ICO 


eoc- 


2 ^- 


dE, 

dx 


lOJ 


2 Hy, 


and another involving Ex, Ey and E^ (H-polarization), 


dH, 

dy 


= -iueoEx 


OH, 

dx 


EqC 


= iujeoE, 




(30) 


(31) 


From Eqs. (27) and (28) we can readily calculate the E-polarized component of the held. Similarly, from 
Eqs. (30) and (31) we will obtain the H-polarized component of the held behind the grating. Once we 
have these expressions, their substitution into Eqs. (22) and (23) leads us to the corresponding EME how 
lines. 

In Figs. andwe can observe a series of sets of EME how lines for different values of the polarization 
angles 0i and 02 (the numerical details of the simulations are provided in the corresponding captions), 
where the wave function just behind each slit, and before they are acted by the polarizers, is given by a 
Gaussian function: 

V’(x', 0+) = -ipiix', 0+) -h 0+), (32) 





















Figure 2: EME flow lines (a) and intensity distribution (b,c) behind two Gaussian slits followed by two 
polarizers and illuminated by monochromatic light with a wavelength A = 943 nm. The parameters 
considered in our simulation have been taken from the experiment data 13 : ai = a 2 = 0.3 mm, 


/ii = —^2 = 2.35 mm, and ai = 02 = l.Scii, and A = 943 nm. The initial polarization is assumed to be 
linear, with A = B, ip = {), 61 = {), and 62 = vr/lO (d). 


where 

V^i(x',0+) = w{x' - pi,ai), (33) 

with i = 1,2 and w{x, a) being the window function defined as 


w(x, a) 


1, X G [— 0 , a] 

0, everywhere else 


(34) 


In Fig. I^a) we notice that the EME flow lines gives rise to a more complete description of interference 
phenomena compared to the description only based on intensity curves (see Eig. [^b)). Specifically, the 
flow lines provide us with an idea on how the electromagnetic field propagates through space, while the 
intensity curves (see Fig.j^b)) only tell us how much intensity is present at each point. In this sense, notice 
how in the far field the maxima and minima displayed by the light intensity correspond, respectively, to 
maximum and minimum values of the density of trajectories. Actually, a slight variation in the relative 
orientation of the axes of the polarizers behind the slits gives rise to non-vanishing minima, as seen 
in right panel of Fig. [^a). Only if the axes of the polarizers are oriented along the same direction, the 
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distribution of EME flow lines is symmetrical with respect to the symmetry axis (see Eig.l^a)). However, 
as their relative orientation increases, as seen in the case displayed in Fig. [^or in Figs. [Sps) and (c), the 
flow-line distribution is not symmetrical and total fringe visibility (vanishing minima) disappears. In this 
regard, the most remarkable case takes place for mutually orthogonal polarizers (see Fig.j^b)), when the 
typical oscillations of intensity disappear. It is in this case when we talk about lack of interference. As 
it can be seen in any of the figures, the distribution of EME flow lines is in compliance with this form of 
intensity -it is uniform, not showing variations (see Fig. j^b)). Apart from the fringe visibility decrease, 
it can also be seen a certain phase shift, downwards for conditions before orthogonality (see Fig. §>)) 
and upwards after it (see Fig. ic)). 

In order to check the theoretically drawn trajectories presented at Figs. and we propose to 
experimentally determine the average photon paths by adding differently oriented polarizers behind the 
slits to the experimental setup reported in 13 . This is a generalization of the previous proposal by 


Davidovic et al. 47 to add orthogonal polarizers in the setup. 


We consider that the proposed experiment could contribute to settle the controversy about the influ¬ 
ence of a potential observer on the form of the interference pattern of two beams with different polar¬ 
izations. The standard interpretation given to the disappearance of interference after inserting mutually 
orthogonal polarizers after the slits is usually based on the Copenhagen notion of the external observer’s 
knowledge (information) about the photon paths, i.e., the slit traversed by the photon on its way to the 
detection screen. Proponents of the principle of complementarity affirm that information on the path 
destroys the interference 

Using EME-flow lines, determined from the EM field and the Poynting vector, one arrives at another 
interpretation. One observes that EME-flow lines starting from slit 1 end up on the left-hand side of the 
screen, while those starting from slit 2 end up on the right-hand side, both with presence of interference 
and with no interference fringes. However, the distribution of EME flow lines is different in each case. For 
the same orientation of the polarizers, the distribution shows interference fringes (see Fig. j^a)), while 
for orthogonal orientations fringes are absent (see Fig. [^b)), in full agreement with the Arago-Fresnel 
laws [^. Hence, the observer’s information on the slit crossed by individual photons (understood as 
energy parcels) seems to be totally irrelevant regarding existence of interference. What really matters 
is the form of the EME field, which eventually influences both the distribution and the topology of the 
trajectories. For values of the mutual angle between polarizers with the interval the symmetry of the 
interferometer setup is broken. Consequently, the intensity distribution is not symmetrical (see Figs. 2 a) 
and (b)). Accordingly, the distribution of trajectories is not symmetrical either, as seen in Figs. 2 a) 
and[^c). Furthermore, a certain number of trajectories starting at slit 1 (covered by a polarizer that 
transmits more energy than the polarizer allocated behind slit 2) cross the symmetry axis, ending up 
on the right-hand side (see Fig. [^a)). Conversely, a certain number of trajectories starting at slit 2 
may cross the symmetry axis, ending up on the left-hand side, if the polarizer allocated behind this slit 
transmits more energy than the polarizer behind slit 1, as seen in Fig. [^c). 


4 Quantum Particle Trajectories - Flow Lines of Quantum-Mechanical 
Currents 

As it has been pointed out, sound waves are described in the framework of mechanics of continuum, even 
though the medium through which they propagate is composed of discrete objects, such as atoms or 
molecules, for instance. In this sense, the pressure, as a fundamental physical quantity in the description 
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Figure 3: EME flow lines behind two Gaussian slits with the same parameters as in Fig. and polarization 
angle = 0 and 02 = 0 (a), 0i = 0 and 02 = 7r/2 (b), and 0i = 0 and 02 = Ttt/IO (c). 


of sound waves, is a continuous function of space and time. This description is well accepted and 
experimentally verified. Nobody opposes or criticizes such a combination of discreteness and continuity. 

In quantum mechanics, wave function T(r, t) synthesizes continuity (wave properties) and discreteness 
(particle properties) in the quantum realm. But, differently from pressure in classical hydrodynamics, 
the physical meaning of the wave function has been open to debate since the very inception of quantum 
mechanics (and, by extension, also quantum optics) .What is the physical meaning of a wave function 
associated with an individual electron, neutron, atom, molecule or photon? This is a longstanding 
question that has been and is still looking for an answer. 
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It is within this context where we consider very useful taking into account the direct analogy between 
the quantum mechanical current density, 


t)V^'(r,t)] = 

Zim 


(35) 


and the Umov vector for sound waves, described by Eq. (0, on the on hand, and Poynting vector for 
electromagnetic waves, dehned by Eq. (23), on the other hand. 


A more complete or, at least, more robust picture of sound wave interference is obtained by con¬ 
sidering flow lines determined by the Umov vector p9 30 . The same happens in the case of photon 


interference when considering electromagnetic energy flow lines 112,42 determined by the Poynting vec¬ 


tor (23). Analogously, a more complete picture of quantum interference with massive particles arises 


after considering flow lines of quantum mechanical probability density 124,25,27,34,39,42,47,50-52 


These more complete pictures might contribute to the resolutions of dilemmas and paradoxes related to 
the question of the physical meaning of the photon wave function as well as for the quantum mechanical 
wave function of nonzero mass particles. 

Using quantum mechanical current density one obtains an objective interpretation of interference 
phenomena as a process of accumulation of single particle events, as conhrmed in experiments with beams 
of one per one, electron 53 , neutron 54 atom and molecule [^. Intensity curves evaluated by taking 


modulus square of a wave function describe only the hnal distribution, obtained after accumulation of 
many particles. Intensity curves do not explain the distributions of tracks of a small number of quantum 
objects. 


5 Concluding Remarks 

The use of flow lines allows us to get a more complete understanding of wave propagation and interference 
phenomena in classical as well in quantum physics. This follows from the simulations and analyses of the 
flow lines associated with sound waves and electromagnetic field in a typical interference device, namely 
a double slit grating. In this way, larger or smaller values of flow-line densities are directly related to the 
common notions associated with the superposition principle of constructive and destructive interference, 
respectively. 

In our opinion, the combined analysis of the propagation and the evolution of flow lines for sound 
waves will be helpful in the long standing debate about the interpretation of the quantum mechanical wave 
function in spite of the different nature of these two kinds of waves. Due to the analogy between classical 
and quantum interference phenomena, the ability to perform in a simpler fashion experiments with 
classical waves that mimic quantum behaviors, should render important insight on quantum systems. 
Notice that sound waves are described by variables that are continuous functions of space and time, 
although gases and fluids, which constitute the material substrate along which the wave propagates, are 
composed of individual, discrete objects, e.g., atoms and molecules. In the same way, trajectories in 
the quantum wave function do not necessarily need to account for the individual motion of a (quantum) 
particle, but provide us information on how probabilities flow in configuration space, and therefore how 
averaged swarms of identical particles travel throughout such a space. 

Based on the results presented here and the recent experiment performed by Kocsis et al. [^, we 
would also like to propose measuring average photon trajectories behind a two-slit grating covered by 
two polarizers for various mutual angles between their polarization axes. We expect that this experiment 
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should render trajectories in compliance with the EME flow lines displayed in Figs.[^and[^ Experiments 
of this kind should contribute to render some light on the longstanding debate on the influence of a 
potential observer on the form of the interference pattern of two beams with different polarizations. 
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